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Abstract. The ESO Slice Project (ESP) is a galaxy red- 
shift survey we have recently completed as an ESO Key- 
Project over about 23 square degrees, in a region near the 
South Galactic Pole. The survey is nearly complete to the 
limiting magnitude bj = 19.4 and consists of 3342 galax- 
ies with reliable redshift determination. 
The ESP survey is intermediate between shallow, wide an- 
gle samples and very deep, one-dimensional pencil beams: 
spanning a volume of ~ 5 x 10 4 h~ 3 Mpc 3 at the sensitiv- 
ity peak (z ~ 0.1), it provides an accurate determination 
of the "local" luminosity function and the mean galaxy 
density. 

We find that, although a Schechter function (with a = 
-1.22, M* ; = -19.61 +51og/iand <t>* = 0.020 h 3 
Mpc~ 3 ) is an acceptable representation of the luminos- 
ity function over the entire range of magnitudes {Mb, < 
— 12 A +5 log h ), our data suggest the presence of a steep- 
ening of the luminosity function for Mb, > —17 +51og/i . 
Such a steepening at the faint end of the luminosity func- 
tion, well fitted by a power law with slope /3 ~ —1.6, is 
almost completely due to galaxies with emission lines: in 
fact, dividing our galaxies into two samples, i.e. galaxies 
with and without emission lines, we find significant differ- 



ences in their luminosity functions. In particular, galaxies 
with emission lines show a significantly steeper slope and 
a fainter M*. 

The amplitude and the a and M* parameters of our lu- 
minosity function are in good agreement with those of the 
AUTOFIB redshift survey (Ellis et al. 1996). Viceversa, 
our amplitude is significantly higher, by a factor ~ 1.6 
at M ~ M*, than that found for both the Stromlo-APM 
(Loveday et al. 1992) and the Las Campanas (Lin et al. 
1996) redshift surveys. Also the faint end slope of our lu- 
minosity function is significantly steeper than that found 
in these two surveys. 

The galaxy number density for Mbj < —16+5 log h is well 
determined (n = 0.08 ± 0.015 h 3 Mpc~ 3 ). Its estimate 
for Mbj < —12.4 +5 log h is more uncertain, ranging from 
n = 0.28 h 3 Mpc~ 3 , in the case of a fit with a single 
Schechter function, to n — 0.54 h 3 Mpc" 3 , in the case of 
a fit with a Schechter function and a power law. The corre- 
sponding blue luminosity densities in these three cases are 
Plum = (2.0,2.2,2.3) x 10 8 h L Mpc" 3 , respectively. 
Large over- and under- densities are clearly seen in our 
data. In particular, we find evidence for a "local" under- 



density (n ~ 0.5n for D c 



< 140 h- 1 Mpc ) and 
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a significant overdensity (n ~ 2n) at z ~ 0.1. When 
these radial density variations are taken into account, our 
derived luminosity function reproduces very well the ob- 
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served counts for bj < 19.4, including the steeper than 
Euclidean slope for bj < 17.0. 

Key words: Galaxies: distances and redshifts - luminos- 
ity function - density 



1. Introduction 

An unbiased and detailed characterization of the luminos- 
ity function of field galaxies is a basic requirement in many 
extragalactic problems. Although several determinations 
of this function are already available in the literature, the 
debate about the faint end slope and the normalization of 
the luminosity function is still open. These quantities al- 
low a "local" normalization which is crucial for the study 
of galaxy evolution and for the explanation of the faint 
galaxy counts. 

Wide angle, shallow samples such as CfA2 (mz < 15.5, 
Marzke et al. 1994), SSRS2 (ra B (o) < 15.5, da Costa et 
al. 1994, Marzke & da Costa 1997) and Stromlo-APM 
(bj < 17.15, Loveday et al. 1992) surveys, can provide 
good determinations of the shape of the luminosity func- 
tion both globally and for different morphological types, 
but the normalization can in principle be significantly af- 
fected by local fluctuations. 

Very deep samples, such as the AUTOFIB composite sur- 
vey (6j < 24, Ellis et al. 1996) and the CFHT redshift 
survey (17.5 < Iab < 22.5, Lilly et al. 1995), are mainly 
designed to study evolutionary effects in the luminosity 
function. 

The ESO Slice Project (ESP, bj < 19.4, Vcttolani et al. 
1997a) is aimed to fill the gap between shallow, wide an- 
gle surveys and very deep, one-dimensional pencil beams, 
allowing a robust estimate of the faint end slope and nor- 
malization of the luminosity function derived from a large, 
uniform and complete sample of galaxies. The redshift dis- 
tribution of the galaxies of this sample, which peaks at 
z ~ 0.1, is deep enough to allow the sampling of a sta- 
tistically representative distribution of the structures in a 
region of the Universe where the evolutionary effects are 
not expected to be important. 

The recently published Las Campanas Redshift Survey 
(LCRS, r ^ 17.5, Lin et al. 1996), based on a large num- 
ber of galaxies with redshift (N ga i <~ 19000), probes es- 
sentially the same redshift depth as the ESP survey, but 
on a larger solid angle. The main differences between the 
LCRS and our survey, as well as most of the existing sur- 
veys, are the photometric band in which the galaxies have 
been selected (red instead of blue) and the fact that the 
low surface brightness galaxies have been deliberately ex- 
cluded from the spectroscopic sample (Shectman et al. 
1996). These facts can make somewhat difficult a direct 
comparison between the LCRS luminosity function and 
others, including ours. 



The organization of the paper is as follows. In Sect. 2 we 
briefly summarize the characteristics of the ESP galaxy 
redshift survey and in Sect. 3 we deal with the problem of 
the K-correction. In Sect. 4 and 5 we describe the methods 
of derivation and the results for the luminosity function 
and the mean galaxy density, respectively. In Sect. 6 we 
discuss our results and compare them with previous esti- 
mates of the luminosity function. In Sect. 7 we summarize 
our results. 

2. The ESO Slice Project 

The ESO Slice Project (ESP) galaxy redshift survey is de- 
scribed in Vcttolani et al. (1997a, hereafter paper I) and 
all the data for the sample will be published in Vettolani 
et al. (1997b, hereafter paper III): here we summarize only 
the main characteristics of the survey. 
The ESP survey extends over a strip of a x S — 22° x 1°, 
plus a nearby area of 5° x 1°, five degrees west of the main 
strip, in the South Galactic Pole region. The position was 
chosen in order to minimize the galactic absorption effects 
(-60° <, b 11 < -75°). The target objects, with a limiting 
magnitude bj < 19.4, were selected from the Edinburgh- 
Durham Southern Galaxy Catalogue (EDSGC, Heydon- 
Dumbleton et al. 1988). 

The right ascension limits are 22 /l 30 m and 01 h 20 m , at a 
mean declination of — 40°15' (1950). We have covered this 
region with a regular grid of adjacent circular fields, with 
a diameter of 32 arcmin each, corresponding to the field 
of view of the multifiber spectrograph OPTOPUS (Lund 
1986) at the 3.6m ESO telescope. The total solid angle of 
the spectroscopic survey is 23.2 square degrees. In order 
to increase the overall completeness, in the last observing 
run we have used the MEFOS spectrograph (Felenbok et 
al. 1997), which has a larger field of view but a smaller 
number of fibers, to observe some of the objects which ci- 
ther had not been observed in the OPTOPUS observing 
runs or for which only poor quality spectra were available. 
We observed a total of 4044 objects, corresponding to 
~ 90% of the parent photometric sample, which contains 
4487 objects: the objects we observed were selected to be a 
random subset of the total catalogue with respect to both 
magnitude and surface brightness. The total number of 
confirmed galaxies with reliable redshift measurement is 
3342, while 493 objects turned out to be stars and 1 object 
is a quasar at redshift z ~ 1.174. No redshift measurement 
could be obtained for the remaining 208 spectra. All the 
details about data reduction, number of objects and com- 
pleteness of the sample are reported in paper I and paper 
III. 

The volume of the survey is <~ 5 x 10 4 h~ 3 Mpc 3 at 
z ~ 0.1, corresponding to the sensitivity peak of the sur- 
vey, and <~ 1.9 x 10 5 h~ 3 Mpc 3 at z <~ 0.16, corresponding 
to the effective depth of the sample. 
The absolute magnitudes are computed as 

M = m-25-5logD L (z) - K(z) (1) 
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where the luminosity distance Dl is given by the Mat- 
tig (1958) expression (throughout the paper we adopt as 
H = 100 km/s Mpc" 1 and q a = 0.5) and K (z) is the 

K-correction. The galactic absorption is assumed to be < K{z) >= ^ fj(z)Kj(z 
negligible, because the strip is close to the South Galactic » 
Pole (Fong et al. 1987). 

Finally, we note that the median internal velocity error for 
our galaxies (~ 60 km/s) is quite small compared to the 
survey depth, and therefore its effect on the calculation of 
absolute magnitudes can be neglected. 



mix at each redshift. This average K-correction is defined 



(2) 



3. The K correction 

Since our database was selected in the blue-green band, 
K-corrections are needed to compute the luminosity func- 
tion even for the moderate redshifts sampled by our galax- 
ies (z < 0.3). The application of the functional forms of 
the K-correction as a function of redshift (e.g. Shanks et 
al. 1984) obviously requires the knowledge of the morpho- 
logical type of each galaxy, which is not available for most 
of our galaxies. We are currently refining a procedure to 
estimate the K-correction for each galaxy directly from 
the spectra, but this work is still in progress. Therefore, 
for this paper we have adopted the following statistical 
approach. 



Fig. 1. Adopted weighted K-correction as a function of red- 
shift. The solid line is a polynomial fit to the points. 



We have computed an average K-correction, as a function 
of redshift, "weighted" with the expected morphological 



where fi(z) is the fraction of galaxies of the i th morpho- 
logical type, at redshift z, whose K-correction is Ki(z). 
The expected fractions of each morphological type, as 
well as their K-corrections, have been computed us- 
ing the Pozzetti et al. (1996) pure luminosity evolution 
model. This model, based on the galaxy spectral library of 
Bruzual & Chariot (1993), is constrained to match the em- 
pirical K-corrections and the morphological mix derived 
locally, and at our moderate redshift the results are al- 
most independent of the details of the adopted evolution- 
ary model. In Fig. 1 we show our "weighted" K-correction, 
computed in bins of 0.025 in redshift (solid circles), with 
a polynomial fit superimposed on the points (solid line). 
The decrease of K(z) for z > 0.2 is essentially due to 
the fact that the expected fraction of early-type galaxies, 
which have a larger K-correction, becomes smaller at this 
redshift, while the fraction of bluer galaxies with smaller 
or almost zero K-correction increases. 
Each galaxy at redshift z' is therefore assigned this 
"weighted" K-correction when computing its absolute 
magnitude. In the non-parametric derivation of the lumi- 
nosity function (see Sect. 4.1) we compute for each galaxy 
the maximum redshift within which it would still be in- 
cluded in the sample (z max ). For this purpose the K- 
correction for z' < z < z max is computed by keeping con- 
stant the morphological mix at the value corresponding to 
the measured redshift z'. In this way the K-correction at 
Zmax can be written as: 



<K{z n 



Z 77 



c) 



(3) 



4. The luminosity function 

4-1. The method 

For most samples of field galaxies the luminosity function 
is well represented by a Schechter (1976) form 



4>(L, x, y, z)dL dV — cf>* 



- L ' L 'd 



dV (4) 



where a and L* are parameters referring to the shape of 
the function and (p* contains the information about the 
normalization; these parameters have to be determined 
from the data. 

Many different methods have been used in the past years 
to compute the parameters of the galaxy luminosity func- 
tion. Recently, however, the STY method (Sandage et al. 
1979) has been the most widely used, and it has been 
shown that it is unbiased with respect to density inhomo- 
geneities (see for instance Efstathiou et al. 1988, Bardelli 
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Fig. 2. a) Normalized luminosity function for 3342 ESP galaxies brighter than Mt, r — —VIA +5 log h . The solid circles are 
computed with a modified version of the C-method (error bars represent lcr Poissonian uncertainties), while the fits are obtained 
with the STY method. Dashed line: single Schechter function; solid line: Schechter function and power law. b) The same as 
panel a), but for galaxies with (open squares and dashed line) and without (filled squares and dotted line) emission lines. For 
clarity only the fit with Schechter function and power law is shown. 



et al. 1991). The basic idea of this method is to compute 
the estimator of the quantity ^ , where $ is the integral 

luminosity function. Under the assumption that the shape 
of the luminosity function is not a function of position [i.e. 
4>(L,x,y, z) dL dV — p(x,y, z)dV ip(L)dL], the probabil- 
ity of seeing a galaxy of luminosity Li at redshift Z{ is 



Pi = 



Mb) 



r+oo 
\ 

JL min {zi) 



(5) 



tp(L)dL 



where L m i n [zi) is the minimum luminosity observable at 
redshift z { in a magnitude-limited sample. 
The best parameters a and L* of the luminosity function 
are then determined by maximizing the likelihood function 
£(a, L*), which is the product over all the galaxies of the 
individual probabilities pi. This corresponds to minimize 
the function 



-21n£ 



N 

a 

I i=l 

N 



^lnLi- N(a + l)lnL*- 



N N , 



T* 



(6) 



where T is the incomplete Euler gamma function and N 
is the total number of galaxies in the sample. 
The normalization parameter <f>* is not determined from 
these equations, and will be derived from the mean galaxy 
density in Sect. 5. 

The STY method is parametric, i.e. it assumes a shape 
for the luminosity function: the derived fit can be com- 
pared with the results of a non-parametric method based 
on the same assumptions. We use a modified version of the 
C-method of Lynden-Bell (1971): given a galaxy of lumi- 
nosity Li at redshift Zi its contribution to the luminosity 
function can be estimated as 



i-£>(£k) 



fe=i 



C+(Li) 



(7) 



where C + (Li) is the number of galaxies with luminosity 
L > Li and redshift z < z max (Li). The resulting integral 
function ip(Li) is normalized to unity at the minimum 
luminosity of the sample. The quantity ip(Li) is computed 
for each galaxy and the results are binned in luminosity 
bins and compared with the Schechter function derived 
with the STY method. 
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Sample M min N ga i a 0* ( h 3 Mpc 3 ) (3 M c 

Galaxies in the -12.4 3342 -1.22t°,°? -19.6ltg;^ 0.020 ± 0.004 
total sample 

-12.4 3342 -1.16 -19.57 0.021 -1.57 -16.99 

Galaxies with -12.4 1575 -1.401^0 - 19 - 47 ±o:ii 0.010 ± 0.002 
emission lines 

-12.4 1575 -1.34 -19.42 0.010 -1.70 -16.94 

Galaxies without -12.4 1767 -0.98±^ -19.621° °* 0.011 ±0.002 
emission lines 

-12.4 1767 -0.90 -19.57 0.012 -1.38 -17.22 



4-2. The results 

The normalized Schcchtcr luminosity function derived 
with the STY method is shown in Fig. 2a (dashed line) 
together with the data points obtained with the non- 
parametric C-method. The error bars on these points cor- 
respond to the statistical (i.e. Poissonian) errors. While 
the Schcchtcr function is an excellent representation of 
the C-method data points for Mbj ^ — 16 +51og/i , at 
fainter magnitudes it lies below all the points down to 
Mbj = —12.4 +51og/i . We have therefore modified the 
model function, adopting a Schechter function for L > L c 
and a power law for L < L c . In this case the low luminosity 
part of the luminosity function is described by: 



(8) 



with [3 and L c being two additional free parameters in the 
fitting procedure, while the normalization A* is fixed by 
requiring continuity of the two functions at L = L c . 
The fit with this two-law function (solid line in Fig. 2a) is 
almost indistinguishable from the single Schechter func- 
tion at bright magnitudes, and reproduces very well also 
the faint part of the luminosity function. Statistically, as 
judged from the decrease of the S function, the improve- 
ment of the two-law fit with respect to a single Schcchtcr 
function is significant at about 2a level, and we can there- 
fore conclude that our data suggest the presence of a steep- 
ening of the luminosity function for Mbj ^ —17 +51og/i . 
In our sample there are 134 galaxies with Mbj > —17 
+51og/i and 38 galaxies with Mbj > — 16 +5 log h : the 
sampled volumes are ~ 1.4x 10 4 h~ 3 Mpc 3 for galaxies 
with M bj = -17 +51og/i and - 3.7 x 10 3 /i~ 3 Mpc 3 for 
galaxies with Mbj = — 16 +5 log h . 

The derived best fit parameters are listed in the first two 
lines of Tabic 1. The errors given for the case of a single 
Schechter function are the projections onto the a and M£ 



axes of the la confidence ellipse (see Fig. 3a) for two in- 
teresting parameters (AS = 2.30). 

The errors for the two-law function are not reported in 
the table because, given the intrinsic correlations among 
the four fitted parameters, it is not very meaningful to 
give the projection on the four axes of the global four- 
dimensional error volume. However, we analyzed all the 
possible combinations of pairs of parameters, considering 
in turn each pair as "interesting" parameters (Avni 1976). 
From this analysis we find that the errors in a and M* are 
very similar to, although slightly larger than, those com- 
puted for the case of a single Schechter function, while M c 
is not well constrained by the data, with a la error larger 
than one magnitude. 

Fig. 3b shows the 2 a contours in the (a, j3) plane, which 
is the most relevant one in determining the reality of the 
steepening at low luminosity. Note that no point of the 
dotted line, corresponding to the locus a = (3 (i.e. no 
steepening) is within the allowed region. 
The faint end steepening is almost completely due to 
galaxies with detectable emission lines. In fact, dividing 
the galaxies into two samples, i.e. galaxies with and with- 
out emission lines (1575 and 1767 galaxies respectively), 
we find very significant differences in their luminosity 
functions. Note that, given the typical signal to noise ra- 
tio of our spectra, detection of a line implies an equivalent 
width larger than about 5 A. Fig. 2b shows the normal- 
ized luminosity functions for galaxies with (open squares 
and dashed line) and without emission lines (solid squares 
and dotted line): for clarity only the fit for the two-law 
function (Schechter function and power law at low lumi- 
nosity) is plotted. It is clearly seen from the figure that 
galaxies with emission lines show a significantly steeper 
faint end slope and a slightly fainter M* (see best fit pa- 
rameters in Table 1). The volume density of galaxies with 
emission lines is lower than the volume density of galaxies 
without emission lines at bright magnitudes, but becomes 
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Fig. 3. a) Confidence ellipses at la and 2a levels for the parameters a and M* , in the case of a fit with a single Schechter 
function, for the total sample and for galaxies with and without emission lines, b) Confidence ellipse at 2a level in the (a, j3) 
plane, in the case of a fit with Schechter function and power law for the total sample. The dotted line corresponds to the locus 
a = /3. 



higher at faint magnitudes. The difference between the 
luminosity functions of galaxies with and without emis- 
sion lines was noted quite early in the project (see for in- 
stance Vettolani et al. 1992), when only a fraction of the 
redshifts was available. Fig. 3a, which shows the la and 
2a confidence ellipses of the parameters a and M* (sin- 
gle Schechter function case) for the total sample and for 
galaxies with and without emission lines, clearly demon- 
strates that the difference between the two subsamples is 
highly significant. 

A similar difference in the best fit parameters of galax- 
ies with and without emission lines has been found in the 
LCRS (see Fig. 10 in Lin et al. 1996), although for each 
subsample their best fit slope is significantly flatter than 
the corresponding slope in our survey. 
Finally, in order to check the possible existence of evo- 
lutionary effects also at our moderate redshifts, we have 
divided the ESP galaxies in three comoving distance bins, 
chosen in order to have three samples with a compara- 
ble number of objects (D com < 250 h^ 1 Mpc ; 250 hr 1 
Mpc < D com < 360 h' 1 Mpc ; D com > 360 h' 1 Mpc ). 
The la confidence ellipses of the (a, M*) parameters 
for the three derived luminosity functions are reported in 
Fig. 4. From the overlap of the ellipses we can say that no 
evidence for evolution is seen in our data, in agreement 



with the conclusion of Ellis et al. (1996) that the bulk of 
the evolution sets in beyond z ~ 0.3. 
As a final check, in order to estimate the maximum 
amount of uncertainty induced by our use of a statisti- 
cal K-correction, we have also computed the parameters 
of the luminosity function with the two extreme assump- 
tions that all galaxies are either ellipticals or spirals. Even 
in these cases, the derived parameters are not much differ- 
ent from those listed in Table 1. The maximum variations 
are Aa ~ -0.06 (i.e. steeper slope) and AM* ~ -0.22 
(i.e. brighter M*) when we apply to all galaxies the K- 
correction appropriate for elliptical galaxies. Moreover, 
in both cases the difference between the parameters for 
galaxies with and without emission lines remains highly 
significant. 



4-3. The influence of magnitude errors and redshift incom- 
pleteness. 

The observed luminosity function is a convolution of the 
true luminosity function with the magnitude error distri- 
bution (Efstathiou et al. 1988). Assuming that the dis- 
tribution of the magnitude errors is Gaussian, with dis- 
persion om (independent of m), the observed luminosity 
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survey with 1/n sampling). If this is not the case, it has 
been shown by Zucca et al. (1994) that eqs.(^) and (6) 
have to be modified as: 



Fig. 4. Confidence ellipses at lcr level for the parameters a 
and M* of the luminosity functions in three different distance 
bins. The asterisk refers to the best fit parameters for the total 
sample. 



function </> j, a 

tptrue as 
4>obs{M) 



is related to the true luminosity function 



-t-oo 



trae (M')e-( M '- M ) 2 / 2 ^dM'(9) 



Preliminary analysis of CCD data in our survey area for 
about 80 galaxies in the magnitude range 16.5 < bj < 
19.4 shows a linear relation between bj(EDSGC) and 
7TIb{CCD) over the entire range of magnitude, with a 
dispersion (<tm) of about 0.2 magnitudes around the fit 
(Garilli et al. in preparation). Since the CCD pointings 
cover the entire right ascension range of our survey, this 
<jm includes both statistical errors within single plates 
and possible plate-to-plate zero point variations. With 
the conservative assumption that the observed ctm = 0.2 
is entirely due to errors on EDSGC magnitudes, we have 
computed the parameters for the case of a Schechter func- 
tion using the convolved luminosity function of eq.([j|) in 
the likelihood function. The resulting best fit parameters 
are very similar to those shown in Table 1 , with the slope 
a becoming flatter by ~ 0.05 and M* fainter by ~ 0.10. 
Even more negligible are the effects on the best fit pa- 
rameters of our small redshift incompleteness. Eq.(H|) is 
correct only for a complete, unbiased sample in which all 
galaxies with to < mu m are members of the sample or all 
galaxies with to < mu m have the same probability of be- 
ing members of the sample (as, for example, in a redshift 



P, 



and 



S = -2 



A' 



i>{L)dL 



(10) 



JV 



< Wi In Li — (a + 1) In L* 



if^-f^nP^ + l,^^) 

i=l i=l ^ ~ ' 



(11) 



where Wi is the inverse of the probability that the i th 
galaxy has of being included in the sample. 
For the ESP sample there are two kinds of redshift in- 
completeness: objects which have not been observed and 
objects whose spectra were not useful for redshift determi- 
nation. We have verified that the former incompleteness 
(~ 10%) is consistent with being random in magnitude. 
As such, it does not affect the derived parameters a and 
M* , but has to be taken into account when determining 
the luminosity function normalization (see next section). 
On the other hand, the latter incompleteness (~ 5%) is 
higher for fainter objects and is well described by: 





0.023(m 



16. 2) 



if to < 16.2 
if to > 16.2 



(12) 



Note that the maximum fraction of galaxies for which the 
spectra did not provide a useful z determination is ~ 7% 
for the faintest galaxies of our survey (bj = 19.4), for 
which the corresponding weight Wi = 1/[1 — /(wi;)] is 1.08. 
Applying this correction to ESP galaxies, we find that the 
effect on the best fit parameters are completely negligible, 
i.e. Aa - A/3 0.02, AM* - AM C ~ 0.01. 

5. The mean galaxy density 

5.1. The method 

Given a magnitude limited sample, an unbiased estimator 
for the mean number density of galaxies n is 



F(z)W(z)^-dz 
dz 



(13) 



where W(z) is a weighting function and F(z) is the selec- 
tion function of the sample, defined as 



F(z) = 



max[Li ■ l L rri i n {z)] 



4>{L)dL 



(14) 



cf>{L)dL 
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where L\ is the minimum luminosity considered in the 
estimate of the luminosity function. The selection func- 
tion represents the ratio between the number of galaxies 
detectable at redshift z and the total number of galax- 
ies with L > L\ at the same redshift. Davis & Huchra 
(1982; DH82) have discussed a number of different esti- 
mators for h, corresponding to different choices for the 
weighting function W(z) 1 and have shown that the min- 
imum variance estimator is nj 3 , obtained with W(zi) = 
1/[1 + hJ 3 F(zi)], where J 3 is the second moment of the 
spatial correlation function £(r). Adopting J3 <~ 10,000, 
as found from a fit of the correlation function in various 
surveys (see for instance Stromlo-APM and LCRS), in- 
cluding ours, and using the formula given in the appendix 
of DH82, we find that the expected fractional error on n is 
~ 10%. Note, however, that this error estimate accounts 
only for the part due to galaxy clustering and does not 
take into account the additional uncertainty arising from 
the errors on the parameters of the luminosity function. 
By varying a and M* along their 1 a confidence ellipse, 
we find that the errors on h induced by the uncertainty 
on the parameters of the luminosity function are signifi- 
cantly larger (<~ 20%) than those due to galaxy clustering. 
Therefore, since there is no a priori reason to prefer the 
nj 3 estimator with respect to other possible estimators, 
we have analyzed the behaviour of all the estimators dis- 
cussed by DH82. 



Fig. 5. Integrated average density as a function of the max- 
imum comoving distance for the four DH82 estimators (see 
text). 



Fig. 6. a) Galaxy density as a function of the comoving dis- 
tance: distance shells have been chosen in order to have a num- 
ber of expected object > 150; error bars represent la Poisso- 
nian uncertainties. The solid straight line corresponds to the 
mean value n, dashed lines represent la uncertainties. The 
values refer to galaxies with M b j < —12.4 +5 log h , in the 
case of a fit with a Schechter function, b) Comoving distance 
histogram of ESP galaxies: the solid line is the distribution 
expected for a uniform sample. 



Fig. 5 shows the estimated average density as a function 
of the maximum distance for the four DH82 estimators: 
the values refer to galaxies with M\, } < —12 .4 +5 log h , 
in the case of a fit with a Schechter function. The relative 
behaviour of the curves does not change with other choices 
for M min or for the case of the fit with a Schechter func- 
tion and a power law. The nz estimator produces a very 
stable value for the galaxy density for all comoving dis- 
tances <; 400 h^ 1 Mpc , even for distances larger than the 
maximum distance in the sample. The difference between 
the maximum and minimum values of the 713 estimator for 
D com > 400 h^ 1 Mpc is of the order of 6%, well below the 
20% uncertainty due to the errors on a and M* . All the 
other estimators appear to be much more sensitive to local 
fluctuations and vary significantly with the assumed max- 
imum distance. For example, the density resulting from 
the nj 3 estimator would change by ~ 25% by changing the 
maximum distance from ~ 650 hr x Mpc , within which 
~ 99% of our galaxies are contained, to ~ 780 h^ 1 Mpc , 
which corresponds to the maximum distance in the sam- 
ple. Since an optimal choice for the maximum distance to 
adopt in this analysis is not defined, this result shows an 
intrinsic ambiguity, at least for our sample, connected to 
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the nj 3 estimator. On the basis of this figure we therefore 
decided to adopt the 713 estimator, which corresponds to 
a constant weighting function in cq.([l3]). The same esti- 
mator has been used in recent analyses of other redshift 
surveys, such as the LCRS (Lin et al. 1996) and the SSRS2 
(Marzke & da Costa 1997). The errors on n are computed 
applying the error propagation on eq.([l3|), taking into ac- 
count also the uncertainties arising by varying a and M* 
along their confidence ellipse. 

The derived densities need to be corrected for the red- 
shift incompleteness of our survey: indeed, among the 4487 
objects of the photometric catalogue, we have 208 failed 
spectra with low signal to noise ratio and 443 not ob- 
served objects (see paper I and Sect. 4. 4). We assume that 
the failed spectra correspond to galaxies, because stars, 
being point-like objects, have on average a better signal 
to noise ratio than galaxies, and that the percentage of 
stars in not observed objects spectra is the same as in 
the spectroscopic sample (i.e. ~ 12.2%). After subtract- 
ing the expected stellar contamination, we are left with 
about 600 not observed galaxies, corresponding to ~ 15% 
of the total. Therefore the galaxy counts are multiplied 
by the factor 1.15, before deriving the mean density. The 
use of a simple multiplicative factor implicitely assumes 
that the distance distribution of non observed galaxies is 
the same as that of the whole sample. This assumption is 
supported by the analysis of the magnitude distributions 
discussed in Sect. 4. 3. 

Finally, from the estimated n, it is possible to derive the 
normalization of the luminosity function as 



(15) 



r 



1 



Li 
L* 



in the case of a single Schechter function, and as 

n 



a-0 



-Lc/L* 



Although the densities derived from most of the distance 
shells are within ±2cr from the global mean density of the 
sample, at least three regions have densities which differ 
significantly (^ a factor of two) from the mean density. 
These regions (two underdense regions at D com < 140 
h^ 1 Mpc and D com ~ 230 h" 1 Mpc , and an overdense 
region at D com ~ 290 h^ 1 Mpc ) are clearly visible also in 
Fig. 6b, which shows the observed distance histogram and 
the distribution expected for a uniform sample. The rela- 
tively large density fluctuations seen in our data are clearly 
due to the fact that our survey is a slice with a narrow 
width in one direction. Note, however, that over-densities 
and under-densities of about a factor two are seen also 
in wider angle surveys, such as the SSRS2 (Marzke & da 
Costa 1997), the CfA2 (Marzke et al. 1994) and the LCRS 
(Lin et al. 1996) surveys. 

The results shown in this Fig. 6a refer to the case of a fit 
with a single Schechter function for galaxies brighter than 
M\,j = —12 .4 +5 log ft . Qualitatively, the results for dif- 
ferent limiting absolute magnitudes or for the case of the 
fit with a Schechter function and a power law are similar to 
those shown in Fig. 6a, except for different normalizations. 
The galaxy number density for M^j < —16 +5 log ft, is 
well determined (n = 0.08 ± 0.015 ft 3 Mpc~ 3 ) and is 
essentially independent from the adopted fitting law for 
the luminosity function. Its estimate for Mbj < —12.4 
+5 log ft is more uncertain, ranging from h = 0.28 ± 0.05 



ft 3 Mpc 



in the case of a fit with a single Schechter 



in the case of Schechter function and power law. 
5.2. The results 

The density radial profile of ESP galaxies is shown in 
Fig. 6a as a function of comoving distance. In each distance 
bin the density has been computed using the 713 estimator 
as discussed above; the error bars represent la Poissonian 
uncertainties. The width of the various distance shells (see 
horizontal lines in the figure) has been chosen in order to 
have ~ 150 expected galaxies in each bin. The solid line 
represents the value of the global n derived from the total 
sample (see Fig. 5) and the dashed lines indicate the ±lcr 
uncertainty on this value. The median of the n values in 
the shells is 0.29 ± 0.04 h 3 Mpc -3 , in excellent agree- 
ment with the global value n = 0.28 ± 0.05 h 3 Mpc~ 3 . 



function, to n = 0.54 ± 0.10 h 3 Mpc~ 3 , in the case of a 
fit with a Schechter function and a power law. The cor- 
responding luminosity densities in these three cases are 
Plum = (2.0,2.2,2.3) x 10 s ft L Mpc" 3 , respectively: 
the similarity of these values indicates that the galax- 
ies in the faint end of the luminosity function contribute 
strongly to the number density, but change only slightly 
the global luminosity density. 

From these number densities we have then derived <fi* for 
ous cases: the obtained values are reported in Ta- 
— blej 1, in the case of M min = —12 .4 +51og/i ; the values 
do not change with other choices for M m i n . 



6. Discussion and comparison with previous re- 
sults 

In Table 2 we list the parameters from the most recent 
luminosity function determinations available in the liter- 
ature. From this table it is clear that there are significant 
differences between the various surveys. We first note the 
discrepant value of M* found in the CfA2 survey (Marzke 
et al. 1994) with respect to all the other samples, which 
cannot be entirely explained by the different passband 
used. The a and M* parameters derived from the SSRS2 
survey (Marzke & da Costa 1997) are in good agreement 
with our results. However, the amplitude of the SSRS2 lu- 
minosity function is such that the ratio between the SSRS2 
and the ESP luminosity functions is in the range 0.4 — 0.5 
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Fig. 7. a) Comparison of the ESP luminosity function with previous results, b) Ratio between the ESP luminosity function 
and the others. In this panel, in order to be consistent with the other luminosity functions, we have used for the ESP galaxies 
the fit obtained with a single Schechter function. The dotted straight lines represent an approximate ±2<r uncertainty on the 
ratio derived from the uncertainties on the normalizations of the various surveys. For sake of clarity, the two shallow surveys 
(CfA2 and SSRS2) are not plotted. 



Table 2. Parameters of the luminosity function from various samples 



Sample N ga i mu m M min a M* <f>* ( h 3 Mpc 3 ) notes 



ESP 


3342 


bj = 19.4 


-12.4 


-i 9 9+0.06 
r-^_0.07 


_1Q 61+0.06 


0.020 ± 0.004 




CfA2 


9063 


mz = 15.5 


-16.5 


-1.0 ± 0.2 


-18.8 ±0.3 


0.04 ±0.01 




SSRS2 


3288 


ra-s(o) = 15-5 


-14 


-1 16+ 008 


-19.45 ±0.08 


0.0109 ± 0.0030 




Stromlo-APM 


1658 


bj = 17.15 


-15 


-0.97 ±0.15 


-19.50 ±0.13 


0.0140 ±0.0017 




LCRS 


18678 


r ~ 17.5 


-17.5 


-0.70 ± 0.05 


-20.29 ±0.02 


0.019 ±0.001 


< bj -r > = 1.1 


AUTOFIBa 


588 


bj = 24 


-14 


-1 16+ 005 

i - 1D -0.08 


-19.30±g"i| 


0.0245+™ 


0.02 < z < 0.15 


AUTOFIBb 


665 


bj = 24 


-16 


_1 41+012 


-19.65±S:iS 


0.0148±±g;Sg?S 


0.15 < z < 0.35 



for -20 < M < -14. 

In Fig. 7a we compare the luminosity functions derived 
from the various surveys; in order to better visualize the 
differences, we show in Fig. 7b the ratio between the ESP 
luminosity function and the others. 
From this figure we find that: 

i) The AUTOFIB luminosity functions (Ellis et al. 1996) 
in their first two redshift bins, which cover approximately 
the same distance range as our survey, are in good agree- 
ment with the ESP luminosity function. Both AUTOFIB 
luminosity functions are within la from the ESP luminos- 



ity function at all magnitudes. 

ii) The Stromlo-APM (Loveday et al. 1992) and the LCRS 
(Lin et al. 1996) luminosity functions are significantly dif- 
ferent from the ESP luminosity function in both shape 
(both are flatter than the ESP at low luminosity) and 
amplitude (both have a lower amplitude by a factor ~ 1.6 
at M~ M*). 

While a possible explanation for the difference in ampli- 
tude with respect to the Stromlo-APM luminosity func- 
tion is presented later in this section, we are intrigued by 
the significant difference in amplitude with the LCRS lu- 
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Fig. 8. Observed galaxy counts compared with those expected 
from the luminosity function (solid line); the same but tak- 
ing into account the observed radial density variations (dashed 
line) 

minosity function, which is derived from galaxies with a 
redshift distribution very similar to ours. The fact that the 
LCRS galaxies are selected in the red band can in principle 
explain, at least in part, their flatter luminosity function 
on the basis of the fact that at low luminosity most of the 
galaxies have emission lines (see Sect. 4. 2 above) and are 
presumably bluer than the average. The dependence of the 
local luminosity function on color has been recently dis- 
cussed by Marzke & da Costa (1997), who indeed found a 
significantly steeper slope for blue galaxies with respect to 
the red ones in the SSRS2. This, however, should not af- 
fect significantly the amplitude at M ~ M* . It is true that 
in the LCRS low surface brightness galaxies have been ex- 
plicitly eliminated from the redshift survey, and therefore 
from the computation of the luminosity function, but the 
fraction of galaxies eliminated because of this selection ef- 
fect is estimated to be less than 10% (Lin et al. 1996). An 
other difference between the LCRS and the ESP surveys is 
in the fraction of stars in the spectroscopic sample of can- 
didate galaxies (~ 3% in LCRS and - 12% in ESP). This 
can in principle be due to a better star-galaxy separation 
in the Las Campanas photometric data. Alternatively, it 
is possible that a non negligible number of galaxies, prob- 
ably the most compact ones, have been classified as stars 
in the Las Campanas data and therefore have not been 
observed spectroscopically. Obviously, we have no way to 
check if something like this has really happened in the Las 
Campanas data. 



We have used our derived luminosity function to compute 
the expected galaxy counts in the case of a constant den- 
sity and compared this prediction (solid line in Fig. 8) with 
the data in our photometric catalogue (solid circles with 
error bars in the same figure). Note that in this magni- 
tude range (bj < 20) the counts predicted using the two 
different fits (single Schechter function or Schechter func- 
tion and power law) are almost identical, while they differ 
significantly at magnitudes fainter than bj <~ 25. 
It is clear from the figure that, while the counts predicted 
from the luminosity function are reasonably consistent 
with the data for bj ^17, they are significantly higher 
at brighter magnitude. Since our counts, although with 
large statistical errors because of the small area, are con- 
sistent with the global EDSGC (Heydon-Dumbleton et 
al. 1989) counts and the APM counts are consistent with 
the EDSGC ones (Maddox et al. 1990), this effect is es- 
sentially the same which has led Maddox et al. (1990) to 
suggest rapid and dramatic evolution in the galaxy prop- 
erties for z <S 0.2. Such a strong "local" evolution is not 
seen neither in our analysis of the luminosity function in 
three distance intervals (see Sect. 4. 2) nor in the deeper 
AUTOFIB redshift survey from which very little, if any, 
evolution is seen up to z — 0.35. 

It has been suggested that the low APM counts at bright 
magnitude may be due to a magnitude scale error in the 
Stromlo-APM galaxy survey (Metcalfe et al. 1995). The 
possible existence of such an error for the APM and similar 
catalogues is reinforced by the recent analysis of well cali- 
brated Schmidt plates digitized with the MAMA machine 
by Bertin and Dcnncfcld (1997). We can not exclude a 
similar problem in our data, although the admittedly lim- 
ited CCD photometry we have obtained on bright galaxies 
in this area (Garilli et al. in preparation) does not show 
any strong magnitude error for galaxies in the magnitude 
range 16 — 17. 

At least for our data, however, an alternative or additional 
interpretation is suggested by Fig. 6, which shows signif- 
icant fluctuations around the mean density. The dashed 
line in Fig. 8 shows the expected counts computed by re- 
laxing the assumption of constant density and allowing 
the normalization of the luminosity function to change 
with distance as indicated by the data points in Fig. 6. 
In this way the agreement between the predicted and ob- 
served counts is more than acceptable over the entire mag- 
nitude range. In particular, the deficiency of the observed 
counts at bright magnitudes would be due to the "local" 
(D com < 140 h~ x Mpc ) under-density. In this respect, 
it is interesting to note that, as mentioned above, the am- 
plitude of the ESP luminosity function is about a factor 
of two higher than the SSRS2 one, which has been de- 
rived from galaxies with D com <, 140 hr 1 Mpc over a 
much larger area, which includes the region of our survey. 
Moreover, an under-density over a similar distance range 
is seen also in the South Galactic Cap part of the LCRS 
(see Fig.8a and 8c in Lin et al. 1996). Since the LCRS 
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is based on galaxies extracted from a much larger area 
(80 x 6.5 sq.dcg.), over which our survey region is fully 
contained, we are led to conclude that the local under- 
density seen in our data may have a size of more than 
100 hr 1 Mpc in at least the right ascension direction. If 
such an under-density extends with a similar size also in 
the declination direction, it could contribute significantly, 
possibly in addition to magnitude scale errors, to the low 
APM counts of bright galaxies. 

In this framework, the difference in amplitude between 
our luminosity function and the Loveday et al. (1992) lu- 
minosity function based on APM galaxies with bj < 17.15 
could be at least partly explained by the presence of a real, 
giant under-density extending over a significant fraction 
of their survey. If we compute a luminosity function with 
our data but with the same limiting magnitude as in Love- 
day et al. , the parameters we derive are fully consistent in 
both shape and amplitude with those derived by them. It 
has to be stressed, however, that the possibility for such 
an explanation has been considered by Loveday et al., but 
they did not find any evidence for it in their data (see, for 
example, their Fig. 4). 

It is also interesting to note that, since a strong "local" 
evolution is not easily accomodated in the standard evo- 
lutionary galaxy models, most models for the faint galaxy 
counts have often assumed a high normalization of the 
local luminosity function (see for example Guiderdoni & 
Rocca-Volmerange 1990, Pozzetti et al. 1996). We find 
that the normalization of our luminosity function is in 
good agreement with that usually assumed in these mod- 
els and therefore our result gives an a posteriori confirma- 
tion of the somewhat ad hoc assumption adopted in these 
models. 

7. Conclusions 

In this paper we have derived the luminosity function and 
the mean density from the ESP galaxy redshift survey; 
our main results are the following: 

1) Although a Schechter function is an acceptable rep- 
resentation of the luminosity function over the entire range 
of magnitudes {M^, < —12.4 +h\ogh ), our data suggest 
the presence of a steepening of the luminosity function for 
Mbj > —17 +51og/i . Such a steepening, well fitted by a 
power law with slope /3 ~ —1.6, is in agreement with what 
has been recently found by similar analyses for both field 
galaxies (Marzke et al. 1994) and galaxies in clusters (see 
for instance Driver & Phillipps 1996). 

2) The steepening at the faint end of the luminosity 
function is almost completely due to galaxies with emis- 
sion lines: in fact, dividing our galaxies into two samples, 
i.e. galaxies with and without emission lines, we find sig- 
nificant differences in their luminosity functions. In par- 
ticular, galaxies with emission lines show a significantly 
steeper slope and a fainter M* . The volume density of 
galaxies with emission lines is lower than the volume den- 



sity of galaxies without emission lines at bright magni- 
tudes, but becomes higher at faint magnitudes. 

3) The amplitude and the a and M* parameters of our 
luminosity function arc in good agreement with those of 
the AUTOFIB redshift survey. Viceversa, our amplitude 
is a factor ~ 1.6 higher, at M ~ M*, than that found for 
both the Stromlo-APM and the Las Campanas redshift 
surveys. Also the faint end slope of the luminosity func- 
tion is significantly steeper for the ESP galaxies than that 
found in these two surveys. 

4) We find evidence for a local under-density, ex- 
tending up to a comoving distance ~ 140 h^ 1 Mpc . 
The volume probed by the ESP within such a distance 
is smaller than the volume of a typical void with 50 
h^ 1 Mpc diameter, and in principle this observed nearby 
under-density could be due to the specific direction of the 
survey piercing through a local void. Our data do not al- 
low to characterize this low density region in terms of size 
and shape. When the radial density variations observed 
in our data are taken into account, our derived luminos- 
ity function reproduces very well the observed counts for 
bj < 19.4, including the steeper than Euclidean slope for 
bj < 17.0. If this under-density extends over a much larger 
solid angle than that covered by our survey, it could, at 
least partly, explain the low amplitude of the Stromlo- 
APM luminosity function. 

5) A similar explanation can not justify the significant 
difference in amplitude between the ESP and the LCRS 
luminosity functions, because the two samples cover essen- 
tially the same redshift range. One possibility, which has 
however to be verified, is that a non negligible number of 
galaxies are missing from the original CCD photometric 
catalog of the LCRS. 

Given the large number of galaxies, the high degree of 
completeness, the accurate selection criteria and the good 
photometry of the ESP redshift survey, we can conclude 
that the results of our analysis give the best available es- 
timate of both the normalization and the faint end slope 
of the luminosity function in the local Universe. 

Acknowledgements. This work has been partially supported 
through NATO Grant CRG 920150, EEC Contract ERB- 
CBRX-CT92-0033, CNR Contract 95. 01099. CT02 and by In- 
stitut National des Sciences de l'Univers and Cosmology GDR. 
We warmly thank Lucia Pozzetti for her help with the weighted 
K-corrections discussed in Sect. 3. We thank the referee for use- 
ful comments. 

References 

Avni, Y., 1976, ApJ 210, 642 

Bardelli, S., Zucca, E., Vettolani, G, Zamorani, G, 1991, A&A 
248, 354 

Bertin, E., Dennefeld, M., 1997, A&A 317, 43 
Bruzual, G, Chariot, S., 1993, ApJ 405, 538 
da Costa, L.N., Geller, M.J., Pellegrini, P.S., Latham, D.W., 
Fairall, A. P., Marzke, R.O., Willmer, C.N.A., Huchra, J. P., 



E.Zucca et al.: The ESP galaxy redshift survey: II. The luminosity function and mean galaxy density 



13 



Calderon, J.H., Ramella, M., Kurtz, M.J., 1994, ApJ 424, 
LI 

Davis, M., Huchra, J.P., 1982, ApJ 254, 437 [DH82] 
Driver, S.P., Phillipps, S., 1996, ApJ 469, 529 
Efstathiou, G., Ellis, R.S., Peterson, B.A., 1988, MNRAS 232, 
431 

Ellis, R.S., Colless, M., Broadhurst, T., Heyl, J., Glazebrook, 

K., 1996, MNRAS 280, 235 
Felenbok, P., Guerin, J., Fernandez, A., Cayatte, V., 

Balkowski, C, Kraan-Korteweg, R.C., 1997, Experimental 

Astronomy, in press 
Fong, R., Jones, L.R., Shanks, T., Stevenson, P.R.F., Strong, 

A.W., Dawe, J.A., Murray, J.D., 1987, MNRAS 224, 1059 
Guiderdoni, B., Rocca-Volmerange, B., 1990, A&A 227, 362 
Heydon-Dumbleton, N.H., Collins, C.A., MacGillivray, H.T., 

1988, in Large-Scale Structures in the Universe, ed. W. Seit- 
ter, H.W. Duerbeck and M. Tacke (Springer- Verlag), p. 71 

Heydon-Dumbleton, N.H., Collins, C.A., MacGillivray, H.T., 

1989, MNRAS 238, 379 

Lin, H., Kirshner, R.P., Shectman, S.A., Landy, S.D., Oemler, 

A., Tucker, D.L., Schechter, P.L., 1996, ApJ 464, 60 
Lilly, S.J., Tresse, L., Hammer, F., Crampton, D., Le Fevre, 

O., 1995, ApJ 455, 108 
Loveday, J., Peterson, B.A., Efstathiou, G., Maddox, S.J., 

1992, ApJ 390, 338 
Lund, G., 1986, OPTOPUS - ESO Operating Manual No. 6, 

ESO, Garching bei Miinchen 
Lynden-Bell, D., 1971, MNRAS 155, 95 

Maddox, S.J., Sutherland, W.J., Efstathiou, G.P., Loveday, J., 

Peterson, B.A., 1990, MNRAS 247, lp 
Marzke, R.O., da Costa, L.N., 1997, AJ 113, 185 
Marzke, R.O., Huchra, J.P, Geller, M.J., 1994, ApJ 428, 43 
Mattig, W., 1958, Astron. Nachr. 284, 109 
Metcalfe, N., Fong, R., Shanks, T., 1995, MNRAS 274, 769 
Pozzetti, L., Bruzual, G., Zamorani, G., 1996, MNRAS 281, 

953 

Sandage, A., Tamman, G., Yahil, A., 1979, ApJ 232, 252 
Schechter, P., 1976, ApJ 203, 297 

Shanks, T., Stevenson, P.R.F., Fong, R., MacGillivray, H.T., 

1984, MNRAS 206, 767 
Shectman, S.A., Landy, S.D., Oemler, A., Tucker, D.L., Lin, 

H., Kirshner, R.P, Schechter, P.L., 1996, ApJ 470, 172 
Vettolani, G., et al. 1992, in International Symposium on Ob- 
servational Cosmology, G.Chincarini et al. eds., ASP con- 
ference series vol.51, p. 74 
Vettolani, G., et al. 1997a, A&A in press [paper I] 
Vettolani, G., et al. 1997b, A&A to be submitted [paper III] 
Zucca, E., Pozzetti, L., Zamorani, G., 1994, MNRAS 269, 953 



This article was processed by the author using Springer- Verlag 
LT E X A&A style file L-AA version 3. 



1 



0.8 



0.6 



o 
. i — i 
+-> 

o 

0) 

o 

o 

0) 

-+J 

.2P 0.2 

0) 



0.4 







-0.2 



J I L 



J I I L 



J I I L 



0.1 



0.2 



0.3 



Redshift 








-24 -22 -20 -18 -16 -14 -12 

Absolute Magnitude 







i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



[] 



a 

+-> 

. I — I 

w 

o 

a 

6 -4 

o 



[] 



[] 



□ 



□ 



■ V 



r 



lilt 



[] 



.■ / 
/ 



/ 



/ 

/ 

, 

/ 

/ 
/ 



■ 



/ 

/ 



J I L 



[fi 



J I I I I I I I I I I I I I L 



J I L 



-24 



-22 



-20 -18 -16 

Absolute Magnitude 



■14 



■12 



-0.60 



GALAXIES WITHOUT 



-1.00 - 
-1.20 " 
-1.40 " 




-20.00 -19.80 -19.60 -19.40 -19.20 -19.00 

M* 



-0.60 



~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

First distance bin 



-0.80 



■1.00 



Second distance bin 
Third distance bin 



-1.20 



-1.40 



-1.60 



-20.00 




-19.80 



-19.60 



J I L 



■19.40 



J I L 



■19.20 



■19.00 








-24 -22 -20 -18 -16 -14 -12 

Absolute Magnitude 



0.5 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 







CO 
H 



O 



LCRS 
APM 

AUTOFIB (0.02<z<0.15) 
AUTOFIB (0.15 z 0.35) 



J I I I I I I I I I I I I I I I I I I I I I L 



-24 



-22 



-20 -18 -16 

Absolute Magnitude 



■14 



■12 




Apparent magnitude 



